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ABSTRACT 

We present a close companion search around sixteen known early-L dwarfs using aperture masking 
interferometry with Palomar laser guide star adaptive optics. The use of aperture masking allows 
the detection of close binaries, corresponding to projected physical separations of 0.6-10.0 AU for 
the targets of our survey. This survey achieved median contrast limits of AK ~ 2.3 for separations 
between 1.2 - 4 A/D, and AK ~ 1.4 at §A/L>. 

We present four candidate binaries detected with moderate to high confidence (90-98%). Two have 
projected physical separations less than 1.5 AU. This may indicate that tight-separation binaries con- 
tribute more significantly to the binary fraction than currently assumed, consistent with spectroscopic 
and photometric over luminosity studies. 

Ten targets of this survey have previously been observed with the Hubble Space Telescope as part 
of companion searches. We use the increased resolution of aperture masking to search for close or dim 
companions that would be obscured by full aperture imaging, finding two candidate binaries. 

This survey is the first application of aperture masking with laser guide star adaptive optics at 
Palomar. Several new techniques for the analysis of aperture masking data in the low signal to noise 
regime are explored. 



1. INTRODUCTION 

The mass determinations of stars through binary stud- 
ies have provided numerous mass-luminosity benchmarks 
for the testing and calibration of stellar models. Such 
studies have only recently begun for the regime of brown 
dwarfs. 

The empirical calibration of brown dwarf models is 
generally made more difficult by the added dependency 
on age in the mass-luminosity relationship. For ex- 
ample, an object spectroscopically classified as a late- 
M dwarf may be a young brown dwarf or an old star 
just above the hydrogen burning limit. This broadens 
the range of potential physical properties (mass, age, 
composition) that generate the same observable spec- 
trum. Conversely, photometry can only very generally 
reveal the objects' physical properties. Measurements 
of brown dwarf masses through the tracking of binary 
orbits provide the strongest constraints on stellar mod- 
els, "mass benchmarks" that reduce the degeneracy of 
photometric stu dies even for targets with unknown ages 
dLiu et al"1l2008f). Mass me asurements of brown dwarfs 
by iKonopackv et all (|2010[ ) show systematic discrepan- 
cies between models and measurements; late-M through 
mid-L systems tended to be more massive than models 
predict, while one T dwarf system was less massive than 



its model prediction. This collection of mass benchmarks 
grows even more important as brown dwarf models are 
extended to infer the masses of directly im aged planets, 
such system HR 8799 (jMarois et all 120081 ). 

Binary surveys have also begun to turn up interest- 
ing statistical results that may suggest the brown dwarf 
binary formation mec hanism is different than that for so- 
lar type binaries (see iBurgasser et al.l (|2006t ) for a sum- 
mary of results from low mass surveys, including many 
results presented in this section). Few surveys, however, 
have produced results for very low mass binaries, espe- 
cially those with very tight separations (< 3.0 AU). This 
regime of short period binaries is particularly challenging 
to achieve with ground-based direct imaging. 

While the companion fraction of brown dwarfs is pro- 
posed to be low (« 15%) and peaked within a nar - 
row separation range, 3-10 AU (Burgass er et al.l [2008) . 
little conclusive results are known for separations less 
than 3 AU despite preliminary evidence that many ad- 
ditional conrrjanions_are^kelv_to r eside at very close dis - 
tances (iJeffries fc Maxtecfl ()2005l ): iPinfield et all (2003); 
IChappelle et al.l (12005T il 

Over 90% of known very low m ass binaries have less 
than 20 AU (|Burgasser et al.ll2006l) . Competing theories 
of stellar formation aim to explain the observed com- 
panion statistics of brown dwarfs. As a general trend, 
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stars appear to have a binary fraction that decreases with 
mass. This can partly be explained by the decreased 
binding energy of lower mass primaries, and thus a max- 
imum binar y separation that decreases as a func tion of 
total mass |Reid et all (f2000h : iClose et all (pool ). For 
very low mass stars, the companion fraction peaks near 
3-10 AU, and exhibits a significant (and statistically sig- 
nificant) drop at separations beyond 20 AU separation. 
Slightly more than half of known very l ow mass binaries 
lay w ithin this narrow separation range (fBuraassc r et al.l 
2006). On the near side of this peak, the data collected 
is very likely incomplete, where the necessary resolution 
(300 mas) stretches the limitations of HST/NICMOS and 
ground-telescopes with AO alone. What data has been 
collected suggests direct imaging may have missed com- 
panions at very close separations. 

Spectroscopic, spectral morphology, and Laser Guide 
Star AO surveys suggest that very tight binaries within 
3 AU m ay be as p l entifu l as binaries of moderate sepa- 
ration. iBurgasseH (|2007|) has used spectral features of 
unresolved sources to indicate composite spectra, im- 
plying multiplicity. This technique has suggested nu- 
merous early /mid-L dwarfs with potential mid-T dwarf 
secondari es and systems of e qual-m ass L/T transition 
objects. Ueffries fc MaxteJ (|2005l) used sparse radial 
velocity data-sets of very low mass systems to pre- 
dict an additional 17-30% binaries at separations less 
than 2.6 AU . Pho tometri c overluminosity stud ies by 
iPinfield et all (120031 ) and iChappelle et all (|2005l ) have 
also hinted at surprisingly larger binary fractions (up to 
50%) in the Pleiades and Praesepe, though concerns over 
membership contamination and the influence of variabil- 
ity limit the conclusiveness of the results. In each study, 
with the exception of the Burgasser mid-L/mid-T sys- 
tems, very low mass binaries tend towards equal mass 
pairs (q ~ 1) at close separations, just as is the case at 
moderate separation. 

These preliminary results contrast those of previous, 
observationally complete surveys that focused on mod- 
erate and wide separation binaries. Those surveys pre- 
dict that fewer than 3% of very low mass comp anions 
sit at separations closer than 3 AU (|Allenll2~007[ ). This 
discrepancy speaks to the importance of additional, ob- 
servationally complete surveys searching for binaries at 
close separations. 

Non-Redundant Aperture Masking (NRM) on 5-10m 
class telescopes, combined with LGS AO, allows sub- 
diffraction limit resolution observations at contrasts high 
enough to search for most binaries in this potentially 
fertile, unresolved region. The detection of close brown 
dwarf binaries, with a typical period of 1-2 years, also al- 
lows mass measurements of late-L or T dwarfs, providing 
particularly valuable empirical benchmarks for the study 
of low mass stellar models. To put into perspective the 
dearth of benchmarks, the mass measurements of fifteen 
very low mass systems (including six with L or T dwarf 
compo nents) using LGS AO alone by iKonopackv et all 
(2010) has tripled the number of very low mass systems 
with mass measurements. 

In Section [2] we describe the sixteen field L-dwarf tar- 
gets imaged at Palomar using aperture masking with 
laser guide star adaptive optics and outline the data anal- 
ysis techniques used to determine the binarity of the tar- 
gets. In Section [3l we present the results of our survey, 



which operated in the range of 60-320 mas (1.1-8.4 AU 
@ 18.4 pc, the median distance of our targets). We iden- 
tify four new candidate L dwarf-brown dwarf binaries 
at moderate or high (90-98%) confidence. This survey 
achieved median contrast limits of AK~2.3 between 1.2 
X/D and 4 A/ D, ruling out companions down to approx- 
imately .06 M for old (5 Gyr) systems and .03 M for 
young (1 Gyr) systems. In Section [4] we discuss the 
aperture masking techniques employed in this paper and 
present recommendations for future faint target obser- 
vations. In Section [5j we summarize the results of this 
survey and discuss its implications for future companion 
searches around brown dwarfs. Finally, we discuss in de- 
tail our methods for calculating detection confidence and 
contrast limits in the Appendix. 

2. OBSERVATIONS AND DATA ANALYSIS 
OBSERVATIONS 

We observed our target sample of sixteen field L dwarfs 
in September and October 2008 with the Palomar Hale 
200" telescope (refer to Tabled]). 

Ten of the sixteen targets in this survey have been ob- 
served previously as part of various co mpanion searches 
using the Hubb le Space Telescope (|Reid et al.l [2006; 
IBouv et al.ll200l . These previous observations were ca- 
pable of resolving low contrast or distant (beyond about 
300 mas) companions. Aperture masking complements 
these previous surveys, extending the detection limits 
around these targets to dimmer and closer companions. 

Aperture masking obse rvations were obtain ed using 
the PHARO instrument (jHavward et al.ll200ll ). with a 
9-hole aperture mask installed in the pupil plane of the 
Lyot-stop wheel (Figure [1]). The longest and shortest 
baselines, which set the approximate inner and outer 
working angle, are 3.94m and 0.71m respectively (58 and 
320 mas in K band). We operate to minimize atmo- 
spheric and AO variation during a single image, using 
PHARO in 256 x 256 sub-array mode with a total of 
16 reads (sub-frames) per array reset and 431 ms expo- 
sures. Every read was saved to disk. In post-processing, 
we discard the first three sub-frames of each exposure 
(usually corrupted by detector reset), and combine the 
remaining sub-frames by a Fowler sampling algorithm in 
which later sub-frames are weighed more heavily. Ap- 
proximately 300 images (each with 16 sub-reads) were 
taken in K s for each target, for a total integration time 
of roughly 60-70 minutes per target. 

The Palomar laser guide star adaptive optics system 
(|Roberts et al.1 120081 ) provided the wavefront reference 
for high order AO correction, while nearby (a few arc- 
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Fig. 1. — The aperture mask inserted at the Lyot Stop in the 
PHARO detector. Insertion of the mask at this location is equiva- 
lent to masking the primary mirror. 



minutes) field stars were observed contemporaneous to 
provide tip-tilt correction. 

Aperture masking operates most effectively when ex- 
posure times are as short as possible, but long enough 
to observe fringes over read noise. The optimal expo- 
sure time depends on the brightness of the targets and 
the level of correction provided by the AO system. Poor 
correction favors shorter exposure times, where variation 
of the incoming wavefront quickly degrades the average 
transmission of long baseline frequencies. For targets 
brighter than about tenth magnitude, the read out lim- 
ited exposure time of the PHARO detector, 431 ms for 
the 256 x 256 array, is sufficient to observe long baseline 
fringes. The targets of this survey are approximately 
twelfth magnitude, and initial experimentation showed 
that short exposures did not consistently provide long 
baseline fringes. Longer exposure times (1 minute) faired 
poorly because variations in correction over the exposure 
degraded the average transmission of long baselines be- 
low background levels. We opted to use short exposures 
and to weigh more heavily in post-analysis those obser- 
vations in which long baseline fringes could be seen (see 
additional discussion later in this section). 

Background subtraction is necessary for targets as faint 
as L Dwarfs and background levels were often compara- 
ble to the signal levels. In many instances our obser- 
vations were background limited. To remove the back- 
ground in post-processing, each target was dithered on 
the 256 x 256 sub-array. 

A requirement for obtaining good contrast limits 
around bright targets is the contemporaneous observa- 
tion of calibrator sources: single stars which are nearby 
in the sky and similar in near-infrared magnitudes and 
colors. This calibration is necessary to remove non- 
stochastic phase errors introduced by primary mirror im- 
perfections and other non-equal path length errors. This 
error can be on the order of one to a few degrees, com- 
parable to the measurement scatter of the closure phases 
for bright targets. For brighter targets, the typical ob- 
serving mode is to obtain several observations of the sci- 
ence target, interspersed with observations of calibrator 
stars. However, the lengthy time of acquisition for the 
laser guide star AO system made this method inefficient 
for this survey. Furthermore, the measurement scatter 
for the faint targets of this survey were much larger than 
the expected systematic error. Therefore, we did not use 
calibrator stars. We note that cali brator stars have als o 
not been used for similar reasons in iDupuv et afl ([2009h . 

APERTURE MASKING ANALYSIS AND 
DETECTION LIMITS 

EXTRACTING CLOSURE PHASES FROM RAW 
IMAGES 

The core aperture masking pipeline implemented in 
this paper is similar to that discussed in previous 
work (|Llovd et al.ll2006t iPravdo et aT]|2006t iKraus et al.1 
l2008f >. with additions to handle low signal to noise data 
and calculate confidence intervals and contrast limits. 

Raw images are first dark subtracted and flat-fielded, 
bad pixels are removed, and the data is windowed by 
a super-Gaussian (a function of the form exp(— fcx 4 )). 
This window limits sensitivity to read noise and acts as 
a spatial filter. A per-pixel sky background map is then 
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Fig. 2. — Interferogram and power spectrum generated by the 
aperture mask. (Left) The interferogram image is comprised 
of thirty-six overlapping fringes, one from each pair of holes in 
the aperture mask. (Right) The Fourier transform of the image 
shows the thirty-six (positive and negative) transmitted frequen- 
cies. (Right, inset and overlay) Closure phases are built by adding 
the phases of 'closure triangles': sets of three baseline vectors that 
form a closed triangle. 

constructed from the set of target data and subtracted. 
The background map is generated by masking out the 
target from each image within a set, then, for each pixel, 
using the median value of the pixel flux from those images 
that were not masked. 

The point spread function of the nine hole mask con- 
sists of thirty-six interfering fringes, called the interfero- 
gram. Because the mask is non-redundant, each fringe is 
produced uniquely by the pairing of two holes; the am- 
plitude and phase of this fringe translates directly to the 
complex visibility of the corresponding spatial frequency. 

Fourier-transforming each image reveals seventy-two 
patches of transmitted power we call splodges (thirty-six 
frequencies transmitted, positive and negative) (Figure 
[2]). The complex visibilities are extracted by weighted 
averaging of the central nine pixels of each splodge. Op- 
tical telescope aberrations, AO residuals, and detector 
read-noise contribute noise to the complex visibilities. 
Under the best conditions, visibility amplitudes suffer 
large (> 100%) calibration errors and are not used for 
the analysis in this survey. 

Visibility phase suffers less from these variations, but 
the use of the comple x triple product and closure phase 
(|Lohmann et al.lll983[ ) yields an observable that reduces 
the effect of wavefront-degradations from baseline- length 
independent sources such as low-order AO residuals. 
For an interferometric array (or aperture mask), closure 
phases are built by adding the visibility phases of 'clo- 
sure triangles': sets of three baseline vectors that form a 
closed triangle (see Figure [2]). The set of closure phases 
have lower noise than visibility phases, allowing precise 
photometric measurements despite the loss of photons 
imposed by the mask. 

Thirty-six baselines are present with the 9-hole mask, 
from which 84 closure phases can be constructed. How- 
ever, these closure phases are not all linearly indepen- 
dent, and the 36 baseline phases cannot be uniquely de- 
termined. The phase information cannot be uniquely in- 
verted (by inverse Fourier transform) into an image with- 
out further assumptions ( see, for example, the CLEAN 
algorithm (H ogboml I1974D an d the Maximum Entropy 
Method (|Gull fc Skillindll984Tl L 

As our survey is a search for binaries, the closure phase 
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signal of such a target can be modeled easily. Thus, 
instead of inverting the closure phases to form an image, 
we search for the modeled binary configuration that best 
fits the measured closure phases. 

TYPICAL RESULTS ON BRIGHT TARGETS 

Aperture masking with natural guide star adap- 
tive optics has been employed during numerous near 
infrared surveys on the Palomar and Keck tele- 
scopes. (For mass determinatio n s made through or - 
bit tracking seelLlovd et al.1 (120061 ) : [Ireland et all (pOcf : 
iMartinache et all (120061 12008ft and lDupuv etaLl (J2009), 
and lKraus et al.l (|2008l ) for an extensive survey of Upper 
Scorpius.) Aperture masking has also recently begun us- 
age in conjunction with K eck laser guide star adaptive 
optics (jDupuv et al.ll2009ft l. 

The observation of bright targets (K s < 9), such as 
nearby early-M dwarfs, with an aperture mask enables 
the detection of companions of contrast up to 150:1 
(AK S ~ 5.5) at the formal diffraction limit and 20:1 
(AK S ~ 3.3) at |A/D in at Palomar. 

In this regime, non-stochastic phase errors introduced 
by the optical pipeline dominate closure phase errors, 
as well as background flux, wavefront residuals of the 
adaptive optics system, and achromatic smearing of the 
fringes. Typically, these sources contribute errors on the 
order of one degree after calibration. 

NOISE PROPERTIES OF DIM TARGETS 

For each star image, of which we have approximately 
300 for each star, we extract closure phases. 

Because our targets are faint and exposure times are 
short, detector read-outs contribute significant noise in 
the phase and amplitude of the complex visibilities and 
bispectrum. Amplitudes are particularly susceptible to 
calibration errors. Even during high signal to noise con- 
ditions, amplitudes have been seen to fluctuate by up to 
100%, and are not directly used for fitting to model bina- 
ries. However, closure phase data show a clear improve- 
ment in per-measurement signal to noise for increasing 
amplitude. That is, bispectrum with the largest ampli- 
tude tend to have the highest fidelity closure phases. In 
order to pare off bad data and weigh higher signal to noise 
measurements more heavily, we empirically estimate the 
relationship between amplitude and closure phase fidelity 
(Figure EJ. 

This relationship is estimated by binning the set of clo- 
sure phase data by amplitude and calculating the stan- 
dard deviation of each bin. As already described, as the 
average amplitude within a bin increased, the standard 
deviation within the bin decreased. To first approxima- 
tion, this estimates the relationship between amplitude 
and closure phase fidelity. 

The noisiest bins often show closure phase errors ap- 
proaching 180°. Because the closure phase is inherently 
a measurement of the bispectrum phasor, there is a 360° 
ambiguity in the measurement of closure phase. Fur- 
thermore, even if the underlying noise source is Gaussian 
distributed, the distribution of measured closure phases 
approaches a uniform distribution when the standard de- 
viation of the noise source is larger than about 180°. Di- 
rect calculation of the root mean squared deviation under 
represents the variance of the underlying noise source; 



the calculation of the mean depends on the choice of an- 
gle zero-point. The variation within one bin was at times 
large enough to motivate alternative methods for aver- 
aging bispectrum data. 

We adopt a maximum likelihood method to calculate 
the standard deviations of bins and overall closure phase 
mean. We presume the closure phases in each bin are 
drawn from a wrapped normal distribution^]. The stan- 
dard deviation is varied to maximize the likelihood of 
the data in the bin. The same mean is used for ev- 
ery bin, and the mean which maximizes the likelihood 
of the entire data set is data set's overall mean. This 
allows bins to take arbitrarily large standard deviations; 
a wrapped normal distribution with large standard devi- 
ation converges towards a uniform distribution. For bins 
dominated by read-noise or very low signal to noise, this 
method accurately estimates very large standard devia- 
tions and translate that into very low weighting for the 
bin. The overall errors of closure phase sets ranged be- 
tween 6-15 degrees. In addition, this method of paring 
off bad data typically reduced errors by a factor of two 
over other methods. 

Even after employing this data paring, some sets of 
closure phase data contained so much noise that no reli- 
able signal could be discerned. In this case, the closure 
phase was removed from the set of eighty-four closure 
phases further analyzed. For some targets, up to half of 
the closure phases triangles were removed. In these cir- 
cumstances, the Mt>-coverage of the data drop allowed the 
possibility of model aliasing: i.e., that multiple binary 
configurations yield similar closure phase sets and each 
fit the data equally well. When previous observations 
of the target were available, we used this information to 
rule out unlikely fits. When not, we list all fits to the 
data. 

Finally, non-stochastic errors are typically on the order 
of one to a few degrees. This contribution is much smaller 
than the statistical error, and as such overall best fits of 
our data changed very little whether or not we attempted 
to calibrate out this component. 

MODELING THE BINARY FIT, AND THE 
CALCULATION OF CONFIDENCE AND CONTRAST 
LIMITS 

For each target, we attempted to fit the observed clo- 
sure phases with a three-parameter binary model (sepa- 
ration p, orientation 9, and contrast ratio r > 1). The 
best fitting model is the one which maximized the overall 
likelihood of the data. Errors in the parameters are cal- 
culated from the curvature of the log-likelihood surface 
at its maximum. 

The strength of our fits were determined by comparing 
the increase in log-likelihood, Alog£, between a single 
star fit and a binary star fit for our real data set as com- 
pared to many simulated data sets of single stars. If the 
real data set has a much higher value of A log C than the 

2 The wrapped normal distribution is the probability distribu- 
tion function of the wrapped variable 8 = x mod 2tt, given by 
Pw(8) = ^2p(8 + 2nk), where p is the unwrapped probability of the 
unwrapped variable x. The sum is over integer values of k from 
— oo to co. The wrapped normal distribution, denoted by WN is, 

WN{8) = ^ y] exp[ — ^ 2er"? 7!k ^ ] ' w 'th t ne same summation 

limits. 
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(b) Closure Phase 84 
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Fig. 3. — Estimating pcr-measurement weights for three closure phase data sets for target 2M 2238+4353. The data sets have compar- 
atively high- (left), moderate- (center) and very low- (right) signal to noises. (Top) Plot of bispectrum (closure) phase vs. bispectrum 
amplitude. Note that larger amplitude data have smaller phase spreads, and a clear asymptotic mean can be identified in the high and 
moderate signal to noise cases. (Closure phase 43 contains no discernible signal, and would be removed from further analysis.) Low ampli- 
tude bispectra are swamped by read noise, introducing phase errors which are nearly uniformaly distributed. The solid line estimates the 
relationship between per-measurement standard deviation and bispectrum amplitude. (Middle) Closure phase vs. approximate weighting. 
Note that the higher weighted points have lower per-measurement standard deviation. (Bottom) Resulting p.d.f. of the closure phase. 



simulated data sets, we regard the real data set to be 
indicative of a binary star. (For purely Gaussian noise, 
Alog£ is equivalent to A% 2 .) 

We simulate measurements of single stars with iden- 
tical noise properties and ui>-coverage of the candidate 
binary target. The measured closure phases of a single 
star is the sum of three sources: the intrinsic signal of the 
target, which is zero for a single star; noise fluctuations 
from various sources which are described by the standard 
deviations measured on the target; and a non-stochastic 
systematic error component, which we assume is negligi- 
ble compared to the stochastic noise of these targets. (As 
a check, we also estimated the typical systematic contri- 
bution from the measured signal of eight survey targets 
whose best fits indicated high likelihood for being single 
stars. Including this component to simulate single stars 
had little effect on the overall confidence measurements.) 

From this information, we generate ten thousand mock 
measurements of single stars. To each, we fit the three- 
parameter binary model, record the Alog£, and build a 
distribution of A log C that result from single star obser- 
vations. 

Comparing the value of A log C of the data's fit to 
the simulated distribution yields the probability of false 
alarm: the probability that our apparent binary fit is an 
observation of a single star co-mingled with noise. The 
confidence that our target is binary is one minus this 
false alarm probability. 

To calculate our contrast limits, we first add model 
binary signals to the simulated single star data. These 
mock binary signals span a range of separations, orienta- 
tions, and contrast ratios. We fit each, determine the fit 
confidence, and determine, for a given separation, the 
highest contrast ratio (i.e. dimmest companion) that 
would be detected with 99.5% confidence (false alarm 
probability of .005). These calculations are discussed in 
more detail in the Appendix. 

CALCULATION OF BAYES' FACTORS 

As an alternative to confidence measure presented in 
the previous subsection, our group also applied bayesian 
methods to calculate the Bayes' Factor of each fit, i.e, 
the odds by which our data favors binary models. 

Using Bayesian comparison, the binary hypothesis is 



tested by contrasting two probabilities: that the data 
set would arise from a binary target observation, and 
that the set would arise from a single star observation. 
Expressed mathematically, this is: 



Pr( binary | data ) 
Pr( single | data ) 
Pr( star is binary 



Pr{ data | binary 



Pr{ star is single ) Pr( data | single ) 



fl 



The first term on the right hand side is an attribute of 
the survey population - it is the ratio of the companion 
fraction to one minus the companion fraction - and is 
independent on the data. 

The second term is the Bayes' Factor, representing 
the odds by which the data favors one hypothesis over 
the other. These probabilities are marginalized (and 
integrated) over the binary parameters. Whereas the 
maximum likelihood method searches out the set of pa- 
rameters that maximizes the likelihood of the data, the 
Bayesian approach averages over the parameters. 



Bayes' Factor 



Pr( binary | data ) 
Pr( single | data ) 



(2) 



/ Pr(p\ bin.)Pr(9\ bin.)Pr(r| bin. )£ (data |p, 9, r) 



£(data | single) 
Pr(pjbinary) 



( 3 ) 

The quantities Pr( ( o|binary), Pr(#|binary), 
Pr(rjbinary) refer to distributions of the compan- 
ion separation, orientation, and contrast ratio as they 
are presumed known prior to our observation. These 
distributions for very low mass primaries are themselves 
ongoing topics of debate and limited by observational 
incompleteness, particularly in the separation regime 
of our survey. (For a current review of separation and 
mass ratio distributions derived throu gh observational 
studie s, the reader is invited to view iBurgasser et al.1 
(2006).) lAllenl (|2007j ) quantifies the underlying compan- 
ion distributions from the currently available data. We 
ultimately chose to use blind prior distributions (known 
as Jeffreys' priors). These distributions are uniform 
for separation and log-uniform for contrast ratio. We 
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compare this choice to the Allen priors and discuss its 
implications. 

Our survey focuses on close binaries; our observations 
probe between roughly 1 an d 8 AU for s eventy five per- 
cent of our L dwarf targets. lAllenl ()2007f ) concludes that 
the (physical) separation of companions can be char- 
acterized by a log-normal distribution which peaks at 
7.2±{$ AU with a la width of roughly 11±~| AU. This 
uncertainty in the peak and width contributes notici- 
ble variability of the resulting distribution at the sep- 
arations we consider (see Figure SJ. One characteris- 
tic unifying the span of distributions, however, is that 
companions closer than ~ 2 AU are less likely by up to 
an order of magnitude. The is due in part to describ- 
ing the distribution as a log-normal functional. This 
choice is motivated by the sharp drop in companion frac- 
tion observed outward of about 10 AU, and by assum- 
ing a similar drop shortward of a few AU, where ob- 
servational data is incomplete. As Allen states, this re- 
sult is derived without well-defined searches for compan- 
ion s at close separations, an d th e preliminary results o f 
the lJeffries fc Maxtedl ([20051 ) and lBasri fc Reinersl (j2006[ ) 
surveys potentially indicate the presence of a larger num- 
ber of close binaries. We use a uniform prior to avoid the 
bias of the Allen priors, keeping in mind that companions 
with high Bayes' factor at less than 2 AU could indicate 
observational evidence of close companions yet may also 
be biased towards undue significance. 



- 




Fig. 4. — Proposed log- normal distribution of companion sep- 
aration around L dwarf primaries from Allen (2007). The peak 
and width of the distribution have been constrained by previous 
surveys. The most likely distribution (solid line) and one sigma 
distributions (dashed lines) are shown. Despite the constraints, 
the distribution is noticeably uncertain in the region of separa- 
tions searched by our survey. We opt to use a uniform prior for 
our Baycsian analysis, noting that such a prior may over signify 
companions closer than roughly 2 AU as compared to the Allen 
prior. Similarly, a confirmed detection of a close companion could 
indicated this distribution has been incorrectly described as log- 
normal (see text). 

Observational surveys of very low mass systems 
sho w a tendency towards equal mass b i naries (q ~ 
1) (|Burgasser et al.|[200l iReid et al.|[200l I Allen! 12001 . 
The distribution of mass ratios has been roughly charac- 
terized by a power law, p(q) oc g 7 , with 7 ~ 2 — 4 depend- 
ing on the survey. The large exponent of this distribution 
indicates that low mass ratio (low q) systems are highly 
unlikely (and rare). Transforming this to a distribution 



of broad-band contrast ratios (r, with r > 1) requires 
assumptions about target age distributions, bolometric 
luminosity corrections, and mass-luminosity models (see 
lAllen et al.1 (|2005l [2003) for these assumptions applied to 
field stars). We wish our prior distribution not depend 
so highly on these assumptions and rather rely on a few 
basic assumptions. 

The rapid drop of the L dwarf mass-luminosity rela- 
tion (i.e., halve the mass of the star and its luminosity 
drops by much more) implies that ratios of contrast are 
larger than ratios of mass, and that contrast ratios still 
favor unity (i.e., p(r) cx (1/r) 1 with < 7 < 2 — 4). The 
blind prior for a scale independent quantity like contrast 
ratios is p(r) oc 1/r which, conveni ently, has th e desired 
properties. It is worth noting that lAllenl (|2007l Fig. 14) 
carries out the transformation from mass ratio to con- 
trast ratio, finding a distribution that follows roughly 
p(r) cx l/(r log r) for contrasts down to below 100:1. 

Finally, the same methods can be applied to calcu- 
late posterior distributions for p, 9, and r for each data 
set. For a data set with a single best fit, this distri- 
bution yields a p.d.f. describing the best fitting param- 
eters. The parameter values and errors quoted in this 
paper are those derived from maximum fit likelihood, as 
discussed in the previous subsection, and are not drawn 
from Bayesian posteriors. However, we calculate the pos- 
terior distributions to assure both methods give compa- 
rable results. 

3. SIXTEEN BROWN DWARF TARGETS - FOUR 
CANDIDATE BINARIES 

Aperture masking is most sensitive to companions be- 
tween X/2D and 4 X/D, corresponding to angular separa- 
tions of 60-450 mas in K s at Palomar and physical pro- 
jected separations ranging from 0.6-10 AU for the targets 
in our survey. 

Our achieved detection limits for all sixteen targets are 
summarized in Table [3J Our limits remain relatively flat 
at separations beyond X/D, plateauing near AK~2.3 for 
more than half our targets, and decline to roughly 1.4 
magnitude shortward of X/D (See figure [5]). 

We infer the (companion) stellar properties and mass 
ratios to the corresponding magnitude limits using the 
DUSTY models for target ages of 5 Gyr and 1 Gyr (Table 
[1]). At the formal diffraction limit (about 110 mas in K s ), 
companions with mass ratios of .83 for 5 Gyr systems and 
.55 for 1 Gyr systems would be resolved for 50% of our 
targets at a 99.5% confidence of detection (Fig. [6]). 

Our survey found four candidate binary systems with 
detections at 90-99% confidence and Bayes' Factors fa- 
voring the binary model (Table [3]). We summarize and 
discuss these detections below. 

For some targets in our survey, closure phase mea- 
surements constructed from the longest baselines had 
too much noise to extract a useable signal. The re- 
sulting drop in w^coverage can give rise to aliasing of 
the model fits: i.e., multiple binary configurations fit the 
data equally well. When possible, we used previous ob- 
servations of the target to rule out certain aliased fits; 
when not possible, all model fits are listed. 

2M 0036+1821: A companion at separation 89.5 mas 
and 13.1:1 contrast was detected with 96% confidence 
and a Bayes' Factor of 7.8:1. The data also fits an al- 
ternative (alias) binary configuration (p ~ 243 mas and 
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TABLE 1 

The Sixteen Very Low Mass Survey Targets 







R.A. 


Jjecl. 




Distance 


T 


IT 
11 


K 


5 Gyr 




_ ( 
<1 { 


m s /m p ) b 




1 Gyr 


Name 


(J2000.0) 


/ Tonnn n\ 
(JzUUU.UJ 


Spectral Type 




(pc) 




(mag) 


(mag) 


(mag) 


65-105/105 


-450 m. 


is 65-105/105-450 mas 


zM UUlo+oolb 


00 


15 


44.76 


lie 1 c no c 

+35 16 02.6 


L2 


20.7 


± 


3 


.2 a 


13.88 


12.89 


12 


.26 


0. 


.87 


/ 


0.84 


0.68 


/ 


0.60 


oat nnoz? i 1011S1 

2JVL UUoo+1821 1 ... 


00 


36 


16.17 


+ 18 zl 1U.4 


L3.5 


8.76 


± 





.06 


12.47 


11.59 


11 


.06 





.86 


/ 


0.83 


0.61 


/ 


0.58 


ZM UU4o+10o4 1 ... 


00 


45 


21.43 


1 1 a q A a A a. 
+ 10 o4 44. D 


L3.D 


10 


.9 


± 


2 


.1" 


13.06 


12.06 


11 


.37 





.83 


/ 


0.82 


0.57 


/ 


0.54 


9M nn4fi-un7i ^ 


00 


46 


48.41 


-U07 It; 177 

TU / TO If .1 


M9 


30 


.5 


± 


4. 


.l a 


1Q QQ 


1 Q 1 Q 
TO. TO 


12 


.55 





.86 


/ 


u.oo 


0.80 


/ 


0.73 


2M 0131+3801 


01 


31 


18.38 


+38 01 55.4 


L4 


20 


.9 


± 


4. 


.2 a 


14.68 


13.70 


13 


.05 


0. 


.92 


/ 


0.91 


0.74 


/ 


0.70 


2M 0141+1804 


01 


11 


03.21 


+18 04 50.2 


L4.5 


12 


.6 


± 


2 


.7 


13.88 


13.03 


12 


.50 





.87 


/ 


0.84 


0.66 


/ 


0.63 


2M 0208+2542"2... 


02 


08 


18.33 


+25 42 53.3 


LI 


25. 


.3 


± 


1. 


.7 


13.99 


13.11 


12 


.59 


0. 


.84 


/ 


0.82 


0.62 


/ 


0.56 


2M 0213+4444 s i... 


02 


13 


28.80 


+44 44 45.3 


LI. 5 


18. 


.7 


± 


1 


.1 


13.50 


12.76 


12 


.21 





.83 


/ 


0.82 


0.56 


/ 


0.53 


2M 0230+2704 


02 


30 


15.51 


+27 04 06.1 


L0 


32 


.5 


± 


4. 


.0" 


14.29 


13.48 


12 


.99 





.88 


/ 


0.87 


0.78 


/ 


0.75 


2M 0251- 0352 s i 


02 


51 


14.90 


-03 52 45.9 


L3 


12 


1 


± 


1 


.1 


13.06 


12.25 


11 


.66 





.92 


/ 


0.91 


0.75 


/ 


0.71 


2M 0314+1603"!... 


03 


14 


03.44 


+16 03 05.6 


L0 


14. 


.5 


± 


1 


.8° 


12.53 


11.82 


11 


.24 





.82 


/ 


0.80 


0.60 


/ 


0.54 


2M 0345+2540" 2 ... 


03 


45 


43.16 


+25 40 23.3 


LI 


26. 


.9 


± 





.36 


14.00 


13.21 


12 


.67 


0. 


.83 


/ 


0.81 


0.57 


/ 


0.53 


2M 0355+1133"!... 


03 


55 


23.37 


+11 33 43.7 


L5 


12 


.6 


± 


2. 


.7° 


14.05 


12.53 


11 


.53 


0. 


.91 


/ 


0.90 


0.77 


/ 


0.72 


2M 0500+0330"!... 


05 


00 


21.00 


+03 30 50.1 


L4 


13. 


.1 


± 


2 


.6" 


13.67 


12.68 


12 


.06 


0. 


.91 


/ 


0.89 


0.72 


/ 


0.65 


2M 2036+1051"!... 


20 


36 


03.16 


+10 51 29.5 


L3 


18. 


.1 


± 


3 


2 a 


13.95 


13.02 


12 


.15 


0. 


.87 


/ 


0.85 


0.63 


/ 


0.58 


2M 2238+4353 


22 


38 


07.42 


+43 53 17.9 


LI. 5 


21. 


.8 


± 


1 


.6 


13.84 


13.05 


12 


.52 


0. 


.84 


/ 


0.82 


0.57 


/ 


0.54 



Note. — Coordinates and characteristics of the sixteen very low mass targets observed in this sample. Photometry 
is taken from the 2MASS catalog. Spectral types (spectroscopic) and distances are taken from DwarfArchives.org, 
unless otherwise noted. "Distance measurements derived from J-band photometry and Mj/SpT calibration data of 
ICruz et all (|2003l ) assuming a spectral type uncertainty of ±1 subclass. ^Survey detection limits of Table [2] given in 
terms of secondary-primary mass ratio, assuming a co-eval system (same age and metalicity). M asses ratios are derive d 
from the 5-Gyr (first row) and 1-Gyr (second row), solar- metalicity substellar DUSTY models of lChabrier et al.l {2000), 
using J and K ban d photometry. Sl Target previously observed bv lReid et al.l (|2006T ) . S2 Target previously observed by 
IBouv et all (|2003l) 



TABLE 2 

Survey Contrast Limits (AK) at 99.5% Confidence 



AK' 



Primary 65.0 85.0 105.0 125.0 145.0 165.0 185.0 225.0 265.0 305.0 345.0 385.0 425.0 



2M 


0015+3516 





.92 


1 


.52 


1 


.73 


1 


.88 


2 


.07 


2, 


.25 


2. 


.27 


2 


.18 


2 


.03 


1 


.82 


2 


.06 


2 


.06 


1 


.79 


2M 


0036+1821 


1 


.77 


2 


.30 


2 


.52 


2 


.57 


2 


.63 


2 


.74 


2 


.77 


2 


.79 


2 


.71 


2 


.70 


2 


.62 


2 


.67 


2 


.56 


2M 


0045+1634 


2. 


,06 


2 


.61 


2 


.82 


2 


.87 


2 


.90 


2. 


.96 


3 


.01 


3 


.02 


2 


.94 


2 


.93 


2 


.80 


2 


.86 


2 


.84 


2M 


0046+0715 





.62 


1 


.01 


1 


.16 


1 


.29 


1 


.48 


1 


.58 


1 


.66 


1 


.60 


1 


.38 


1 


.27 


1 


.28 


1 


.35 


1 


.22 


2M 


0131+3801 


0. 


,74 


1 


.26 


1 


.30 


1 


.30 


1 


.35 


1 


.47 


1 


,52 


1 


.55 


1 


.45 


1 


.28 


1 


.25 


1 


.41 


1 


.25 


2M 


0141+1804 


1 


,52 


2 


.13 


2 


.37 


2 


.51 


2 


.55 


2 


.61 


2. 


,65 


2 


.59 


2 


.58 


2 


.58 


2 


.41 


2 


.51 


2 


.42 


2M 


0208+2542 


1 


,29 


1 


.93 


2 


.16 


2 


.28 


2 


.35 


2 


.48 


2 


.52 


2 


.47 


2 


.34 


2 


.28 


2 


.29 


2 


.32 


2 


.25 


2M 


0213+4444 


1 


,84 


2 


.40 


2 


.59 


2 


.64 


2 


.72 


2. 


.77 


2. 


,79 


2 


.81 


2 


.76 


2 


.73 


2 


.59 


2 


.69 


2 


.58 


2M 


0230+2704 





,72 


1 


.11 


1 


.20 


1 


.18 


1 


.23 


1 


.27 


1 


.30 


1 


.32 


1 


.27 


1 


.09 


1 


.03 


1 


.26 


1 


.08 


2M 


0251-0352 





,69 


1 


.07 


1 


.24 


1 


.26 


1 


.32 


1 


,39 


1 


.48 


1 


.38 


1 


.36 


1 


.28 


1 


.24 


1 


.37 


1 


.34 


2M 


0314+1603 


1 


,52 


2 


.08 


2 


.32 


2 


.47 


2 


.54 


2 


.60 


2 


.65 


2 


.61 


2 


.54 


2 


.50 


2 


.51 


2 


.52 


2 


.39 


2M 


0345+2540 


1 


,75 


2 


.28 


2 


.51 


2 


.56 


2 


.61 


2. 


.70 


2. 


.76 


2 


.75 


2 


.61 


2 


.57 


2 


.55 


2 


.58 


2 


.51 


2M 


0355+1133 





.69 


1 


.15 


1 


.21 


1 


.05 


1 


.04 


1 


.22 


1 


.27 


1 


.25 


1 


.29 


1 


.32 


1 


.24 


1 


.27 


1 


.12 


2M 


0500+0330 





,79 


1 


.35 


1 


.53 


1 


.64 


1 


.81 


1 


,97 


2. 


.02 


1 


.99 


1 


.80 


1 


.57 


1 


.80 


1 


.83 


1 


.57 


2M 


2036+1051 


1 


.30 


1 


.90 


2 


.10 


2 


.26 


2 


.31 


2 


.41 


2 


.50 


2 


.40 


2 


.31 


2 


.24 


2 


.26 


2 


.30 


2 


.18 


2M 


2238+4353 


1 


,77 


2 


.29 


2 


.51 


2 


.53 


2 


.57 


2. 


,63 


2. 


.71 


2 


.69 


2 


.57 


2 


.54 


2 


.52 


2 


.56 


2 


.50 



Note. — Detection contrast limits around primaries: "Primary-Secondary separations are given in units of mas, 
and the corresponding detection limits are in AK magnitudes. 
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TABLE 3 
Model Fits to Candidate Binaries 





J. Date 


Separation 


Az. Ang. 


Contrast 


Bayes 




Separation 


Primary 


(+245000) 


(mas) 


(deg) 


Ratio 


Factor 


Conf. 


(AU) 


2M 0036+1821 


4731 


89.5 ± 11.4 


114.1 ± 5.5 


13.14 ± 3.14 


7.9 


96% 


0.78 ± 0.10 


2M 0345+2540 


4731 


217.4 ± 9.1 
352.7 ± 10.5 


258.8 ± 2.8 
87.6 ± 2.0 


26.44 ± 4.22 
30.79 ± 9.08 


7.6 


98% 
96% 


5.85 ± 0.26 
9.49 ± 0.31 






128.2 ± 10.3 


209.9 ± 5.3 


17.76 ± 4.25 




97% 


2.79 ± 0.30 


2M 2238+4353 


4732 


228.5 ± 9.1 


251.8 ± 3.5 


23.79 ± 5.92 


7.1 


95% 


4.98 ± 0.42 






395.5 ± 9.7 


19.5 ± 1.2 


17.63 ± 4.22 




97% 


8.62 ± 0.66 


2M 0355+1133 


4757 


82.5 ± 13.0 


276.2 ± 4.1 


2.10 ± 0.40 


6.3 


90% 


1.03 ± 0.27 





Fig. 5. — Contrast limits at 99.5% detection as a function of primary-companion separation: (left) The primary-secondary magnitude 
difference in K s detectable at 99.5% confidence, (right) The same detection limits in terms of the absolute magnitude of the companion. 
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Fig. 6. — Companion mass and mass ratio limits at 99.5% detection as a function of primary-companion separation: (top left) The 
primary-companion mass ratio detectable at 99.5% confidence. Dashed lines are for systems aged 5 Gyr; Dot-dashed lines are systems ages 
1 Gyr. (top right) The same data in terms of companion mass, (middle/bottom left) As a function of separation and companion mass, this 
plot reveals the percentage of 5 Gyr (middle) and 1 Gyr (bottom) companions detectable at 99.5% given the data quality of the survey. 
Binaries in the white area would have been detected for 100% of the survey targets, followed by contour bands of 95%, 90%, 75%, 50%, 
25%, and 10%. At the diffraction limit (110 mas), companions of mass 0.65 Mq would be resolved for 50% of our targets, (middle/bottom 
right) The same plot in terms of mass ratio. Diffraction limit sensitivity: 5 Gyr companions of mass 0.65 Mq (.038 Mq for 1 Gyr) would 
be resolved for 50% of our targets. Equivalently, our survey reached mass ratios of .83 (5 Gyr) and .55 (1 Gyr) for 50% of our targets at 
the diffraction limit. 

25:1 contrast) with 96% confidence that w e rule out by 
a previous observations. iReid et all ([2006) observed this 
target in November 2005 with the NICMOS imager on 
the Hubble Space Telescope in the F170M and F110W 
bands. At or near this separation, this alternative config- 
uration would have likely been detectable in the F110W 
bands. 

2M 0355+1133: A companion at separation 82.5 mas 
and 2.1:1 contrast was d etected with 90% confidence and 
a Bayes' Factor of 6.3:1. IReid et all (|2006l ) also observed 
this target in the F110W band and found no compan- 
ion. As a proxy for the F110W bandpass, we estimate a 
J band con trast of 2.5:1 using the J-K color- magnitude 
relations of iDahn et all (|2002j ). Their program achieved 
a contrast limit of 2.5:1 beyond approximately 100 mas 
in F110W, suggesting that this candidate binary sat at 
the edge of their detection limits. 

2M 2238+4353: Thirty-five percent of the closure 
phase triangles showed very high noise and were removed 
from analysis. As a result, aliasing of the signal was par- 
ticularly problematic. Three distinct binary configura- 
tions were detected at 95-97% confidence. These range 
in separations between 100 and 400 mas and contrasts 
between 17:1 and 28:1. 

2M 0345+2540: Like 2M 2238+4353, a large per- 
centage of the closure phases were removed from anal- 



ysis. Two distinct configurations, both with contrasts 
—28:1 (AK ~ 3.5 ) were determined with high confi- 
dence. IBouv et al.l (|2003T ) observed this target with the 
Wide Field Planetary Camera 2 (WFPC2) on the Hub- 
ble Space Telescope in March 2001, but we estimate these 
companions to be below their detection limits. Their sur- 
vey reached background limitations at contrasts between 
AM -3-5 in the F814W band. Using the I band as 
a proxy for F814W, we estimate the companion of 2M 
0345+2540 to have a contrast of AJ >5 and to have been 
undetectable in the Bouy survey. 

2M0213+4444: We observed target 2M0213+4444 
three times over two nights in September 2008 (two sets 
in K s , one in H) and once one month later (in K s ). Two 
data sets from September were of poor quality and were 
not used for analysis. The remaining set from September 
found one binary fit (p — 81 mas, 9 ~ 290°, 5.2:1 con- 
trast in K s ) at 89% confidence. The target was observed 
again in K s in October under poor seeing and much of 
the data was unusable. This data set could not be fit well 
by the September results, and implied a different config- 
uration with 90% confidence (p — 234 mas, 9 — 135°, 
11:1 contrast in K s ). Given the low confidences of fits 
and the unreproducibility of these results, we conclude 
that this target is unlikely to be binary. 
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4. DISCUSSION: APERTURE MASKING OF FAINT 
TARGETS 

The use of non-redundant masking removes many 
types of spatial perturbations to the incoming wavefront. 
During high signal to noise observations, when read and 
background noise are minimal, the largest contributor 
to measurement noise is the temporal and spatial atmo- 
spheric fluctuations of the wavefront, even after adap- 
tive optics correction. Short exposure times, roughly less 
than the coherence time of atmosphere, freezes the tip- 
tilt and low-order perturbations to the wavefront, which 
can be removed by combining fringe phases into closure 
phases. This advantage is lost when exposures extend 
over multiple coherence times. For this reason, aperture 
masking flourishes with short exposures. 

Behind laser guide star adaptive optics systems, al- 
though the structure of the corrected wavefront may be 
different, the functionality of aperture masking is the 
same. However, targets requiring laser guide star AO 
tend to be fainter, and thus require either longer ex- 
posure times (permitting sufficient correction) or tech- 
niques for dealing directly with noise from read outs and 
background flux. 

This survey opted for maintaining short exposure 
times. The signal to noise of fringe amplitudes decline 
rapidly for longer baselines, as the transmission function 
for these baselines is lower and turbulence variations are 
larger. Just as, for instance, Stehl ratio depends on the 
variance of the incoming wavefront, so does the fringe 
amplitude, also dropping as exp(— o\ aseline ). For faint 
targets, long baselines fringes often linger undetectable 
below the background and read noise, making difficult 
measurements of long baseline phases. 

The capture of a large number of short exposure im- 
ages allows us to select out the best fringe measurements, 
during the serenditous moments of very good correction 
or still atmospheres, and discard those dwarfed by read 
noise. This technique, analogous to lucky imaging, effec- 
tively selects high signal to noise measurements of closure 
phase. In most cases, these lucky closure phases were 
sufficient to obtain measurements of the target closure 
phase, even at long baselines. 

We contrast this method to two measurements of tar- 
gets observed with long (1 minute) exposures. These 
exposures did often have long baseline fringes detectable 
at or just above background. But this method resulted 
in poor measurements of the target closure phase, even 
at shorter baselines. The multiple-coherence time expo- 
sures means that low order perturbations are not effec- 
tively removed by closure phases, resulting in large phase 
errors, and the fewer overall number of data points re- 
moves the statistical advantage. The measured closure 
phase is not a good measurement of the true target phase. 

Long exposures, with adequate correction, do allow 
longer baseline fringes to grow in amplitude above the 
read noise or background limit. Exposures for aperture 
masking should be limited to the effective coherence time 
of the adaptive optics system - the interval over which 
the phase variance of the longest baselines reaches about 
one radian. 

The quality of measurements from both sets of expo- 
sure data suggests a slight modification of technique for 
the next application of aperture masking with laser guide 



star adaptive optics. The higher noise content of the one- 
minute exposures suggests that these exposures are too 
long for the level of correction obtained in this survey. 
The short exposure method fared much better, but a 
large percentage of images failed to observe fringe ampli- 
tudes above read noise. This suggests that slightly longer 
exposures would have benefited the observations. It is 
worthwhile to note that the low transmission of the long 
baselines, even at Strchl ratios of 15% typically reached 
in this survey, indicates that direct imaging would not 
have been able to obtain X/D resolution. 

5. CONCLUSION 

We present the results of a close companion search 
around nearby L dwarfs using aperture masking inter- 
ferometry and Palomar laser guide star adaptive optics. 
The combination of these techniques yielded typical de- 
tection limits of AK S = 1.5-2.5 between 1-4X/D and lim- 
its of AK S = 1.0-1.7 at 0.6 X/D. Our survey revealed four 
candidate binaries with moderate to high confidence (90- 
99%) and favorable Bayes' Factors. 

Ten of the targets have previously been observed 
with the Hubble Space Telescope as part of companion 
searches. As such, we did not expect to find bright or dis- 
tant companions around these targets which would have 
been identified in the previous surveys. But as demon- 
strated in this paper, the detection profile of aperture 
masking is capable of revealing close or dim compan- 
ions which are obscured by the point spread function of 
full aperture imaging. Aperture masking demonstrates 
an increase in formal resolution and detectable contrast 
at close separations over laser guide star adaptive optics 
alone. 

Our survey indicated two previously observed targets 
as candidate binaries. Our survey indicated one com- 
panion around 2M 0355+1133 within the formal diffrac- 
tion limit of the HST and one companion around 2M 
0345+2540 below the background detection threshold of 
the previous survey. Two other targets, 2M 0345+2540 
and 2M 2238+4354, also indicated the presence of com- 
panions, both with contrast ratios greater than 15:1. 

Aperture masking is most sensitive to companions be- 
tween X/2D and 4A/D, corresponding to angular sep- 
arations of 60-450 mas in K s at Palomar and physical 
projected separations ranging from 0.6-10 AU for the tar- 
gets in our survey. Two candidate binaries presented in 
this paper have projected separations less than 1.5 AU. 
The results suggest a favorable target set for future com- 
panion searches. Their candidacy is consistent with the 
conjecture that tight binaries are underrepresented in the 
current tally of low mass binaries. Spectroscopic surveys, 
which focus on separations within 3 AU, are necessary to 
conclusively answer this question. Extending the use of 
aperture masking with laser guide star AO is a reward- 
ing approach for detecting companions within this range, 
and facilitating the measurements of their masses. 
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APPENDIX 

BINARY DETECTION CONFIDENCE AND CONTRAST LIMITS 

In this section, we provide the basic formalism for fitting closure phase measurements to binary models. This 
formalism also drives the Monte Carlo simulation which determines the strength of these fits. 

Each pair of holes in the Palomar non-redundant nine-hole aperture mask is designed to transmit one unique 
Fourier spatial frequency, for a total of thirty-six frequencies transmitted by the mask. Each image produced by the 
mask consists of thirty-six overlapping fringes which, when Fourier transformed yields the amplitude and phase of 
each transmitted frequency. Combining these phases by adding specific triplets to form closure phases produces an 
observable that is more robust to many forms wavefront noise. Eighty-four closure phases are extracted from each 
image, then averaged over the set of images. Finally, these averages are compared to model closure phases of various 
binary configurations to determine the likelihood that the target is binary. 

DETERMINATION OF BEST FIT 

The measured closure phase signal of each image is the composite of three sources: the intrinsic signal of the target, 
which is zero for a single star and non-zero for a binary; a non-stochastic systematic error component, which may vary 
from target to target, but not during the observation of a single target; and stochastic noise from various sources such 
as atmospheric turbulence, read noise, etc. We denote the intrinsic signal by ^ binary, the systematic component by 
(3 sy stem(t)-> an d the stochastic noise by £ n oise(t,i)- That is, the measurement of closure phase k, extracted from image 
i, observed during the target acquisition set t, &k,t,i, is: 

^k,t,i = ^k,binary + Pk,system(t) + £k,noise(i,t) ■ (Al) 

The measured closure phase from a given set of images that will be fit to model binaries is 

,binary sy stem(i) i 

(A2) 

where the sum is taken over the set of images, and the noise properties of the mean is described the by distribution of 
£,k,noise(i) convolved over the set of images. The systematic component Pk, sy stem(t) ma y change from one acquisition 
to another (affected by telescope movement or AO system performance, etc.) but is assumed constant during the 
observation of a single target. 

We typically use the measurement of calibrator (single) stars, with zero intrinsic signal (i.e., ^k,bmary = 0), to 
estimate the underlying distribution of systematic noise in the optical system. The typical observing mode is to obtain 
several observations of the science target, interspersed with observations of calibrator stars. Although the systematic 
component cannot be determined exactly because it is itself a random variable, we can compile a composite distribution 
of Pk, system from calibrator observations. Subtracting the systematic component from the measured closure phases 
(or, specifically, convolving the two distributions), leaves remaining the intrinsic signal of the target, ^k,binary 

(This calibration step is important for obtaining high contrasts during high signal to noise observations, when the 
contribution from systematic noise is on the order of the stochastic noise. Calibration is less effective when the typical 
stochastic noise is much larger the systematic component.) 

We wish find the three-parameter model binary (separation p, orientation 8, and contrast ratio r) which best fits 
the calibrated signal, ^kfiinary Approached as a maximum-likelihood problem, we calculate the probability that our 
measured data would result from a noiseless, modeled binary, 1 i m (p, 6, r), observed through the noise in the system: 

,£>mary}) °^ J^J / £*k(^k,t — ^ m Pk, system) p{Pk, system) ^Pk, system (-^-3) 
k J 

The integral is due to the convolution of the systematic distribution with the likelihood function, which itself follows 
the distribution of $,k,noise- 

The best-fitting model is that which maximizes the above probability, which we determine by a combination of 
gradient search and visual inspection. The parameter errors are calculated from the curvature of the log-probability 
surface at the maximum. Calculating the confidence of this fit, i.e. that this model represents the true target 
configuration, is detailed in the next subsection. 

Assuming the underlying distributions are wrapped-normal or Gaussian, the convolution above reduces to a single 
wrapped-normal or Gaussian distribution with mean < ^k,t > — < Pk, system > and variance a% k t + a^ k systcm ■ The 
maximum probability problem reduces to one of minimizing \ 2 ■ 
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BINARY DETECTION CONFIDENCE 

Our null hypothesis, which we wish to test against the binary fit, is that the observed target is a single star, with 
intrinsic binary phase zero. Following Eq (|A3j) . the probability of the null model is 

Cnull = PnuUi^m = 0|{#k,5i„ O ry}) QC JJ / p(* fc)t = /3 k \f3 k ) P {/3 k ) d(3 k . (A4) 

k J 

A natural goodness-of-test statistics for comparison is the ratio of the data likelihood to that of the null likelihood, 
or the log of this ratio: 

l°g( T° e ) = log(£ mo d e z) - log(£ n „«) = Alog£. (A5) 

Systematic and stochastic noise may at times conspire to mimic a binary signal, as expressed by a higher probability 
of a binary model, even though the target is a single star. This is a false alarm event. We, therefore, classify the target 
goodness-of-fit statistic as statistically significant only if its value is large compared results of noisy observations of 
single stars. 

We simulate ten thousand measurements of single stars with identical (it, i;)-coverage and noise properties of the 
candidate binary target. The intrinsic phase of a single star, ^binary, is zero. For one measurement, the contribution 
due to statistical noise is drawn from the measured distribution of £k,noiset which typically can be approximated 
by a wrapped normal or Gaussian distribution with mean zero and measured standard deviation. The systematic 
contribution, if included, is drawn from a distribution /3 ktS y Stem , compiled from observations of calibrator (single) 
stars. 

We then apply the same approach used for the target star to each simulated single star. We fit the three-parameter 
binary model to each simulated single star, record the A log C, and build a distribution of A log(£) that results from 
fitting single stars. The probability that our target fit is statistically significant, then, is the probability that the 
goodness-of-fit of a single star is less than the target data's goodness-of-fit: 

V false alarm (* ro )=p( (Alog£) best 

fit to data 

< (Alog£)/ij s to single stars ) (A6) 

and 

detection confidence = 1 -p talse alarm- (A7) 

We consider the target data to reveal a definitive binary detection if the best-fitting model produces a detection 
confidence greater than 99.5% (false alarm probability less than 0.5%). Note that this empirical method is more 
conservative than comparing the values of A log C, which reduces to A\ 2 for Gaussian noise, to a distribution of \ 2 
variables with three degrees of freedom (Fig. [7]). 




i1l,log firobQln it ,-) dll(?g — probolil'ityj 



Fig. 7. — Determination of Binary Confidence for Target 2M 0036+1806. The goodness-of-fit statistic, Alog£ = 6.55 (vertical line), 
is compared to a distribution generated from fits to simulated single stars, resulting in a fit confidence of 96% (empirical). Notice that 
comparing this value to a x 2 distribution with three degrees of freedom (Gaussian) results in a much higher confidence of fit. 
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CALCULATION OF CONTRAST LIMITS 

Whether the target observed is identified as a single star or binary, we are able to quantitatively state the highest 
contrast (dimmest) companion that our technique would have been capable of identifying with high confidence (99.5%) 
as a function of separation. This is, in essence, a statement on the noise characteristics of the data and the w^coverage 
of our mask. 

This amounts to asking the following question: Given simulated binary observations (separation p, orientation 0, 
and contrast ratio r), at what contrast does our detection confidence drop below 99.5%? 

Simulated binary data is the composite of the same noise contributions to single star data plus an intrinsic signal 
due to the presence of a companion. That is, the nth simulated binary data is: 

*k,binary = ^single + * k, model {p, 0, r) (A8) 

For binary models of a range of separations, orientations, and contrast ratios, we generate ten thousand mock binary 
signals of each by adding the intrinsic binary signal to the mock noise simulations described in the previous subsection. 
We fit each, record the fit confidence, and determine the average confidence that that binary can be detected under 
conditions identical to the target observation. For each separation, averaged over orientations, we determine the 
highest contrast ratio (i.e. dimmest companion) that would be detected with 99.5% confidence (false alarm probability 
of .005). 

Ideally, we would determine the confidence of each mock binary by search for its best fit, recording its A log C, and 
comparing it to the false alarm distribution of the previous subsection. In practice, using a fitting routine to fit each 
of these simulated binaries is computationally slow. 

Instead we approximate this process by modifying the false alarm distribution. We make the assumption that the 
inserted binary model yields the best fit. Because we effectively restrict the fitting search to the range of separation, 
orientations, and contrast ratios used to generate the mock binaries, we apply the same restriction to the fitting 
search that generates the false alarm distribution. We then use this modified false alarm distribution to determine the 
confidence of the mock binary fits. In practice, this approximation produces contrast limits slightly more conservative 
than full fitting by about 5-10%. 
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